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ABSTRACT 

The lecture is devoted to the comparison of a few models of cosmic recombination kinetics with 
recent CMB anisotropy data and to corresponding predictions for the upcoming Planck mission. 
The influence of additional sources of ionized photons at the epoch of hydrogen recombination 
may restrict our possibility to estimate the density of baryonic fraction of matter when using 
Planck data. This new source of degeneracy of the cosmological parameters which can be 
estimated from the CMB anisotropy data can mimic the distortion of the CMB power spectrum 
for different models of baryonic matter. I would like to point out that realistic values of the 
cosmological parameters can not be obtained without Planck polarization data. 
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1. Introduction 



In my lecture, I will discuss the physics of cosmological hydrogen recombination and properties of the 
last scattering surface (at redshift z ~ 10 3 ), which play a crucial role for the Cosmic Microwave Background 
(CMB) anisotropy and polarization formation. A detailed understanding of the ionization history of the 
cosmic plasma is very important in the context of recent CMB experiments, such as the BOOMERANG (de 
Bcrnardis et al. 2000), MAXIMA-1 (Hanany et al. 2000), DASI (Halvcrson et al. 2001), CBI (Mason et al. 
2002), VSA (Watson et al. 2002) and especially, the new DASI polarization data (Kovac et al,2002 ). I will 
show how these experiments, including the upcoming MAP and Planck data, would be fundamental for 
our understanding of the most general properties of the structure and evolution of the Universe, the history 
of the galaxies and the large-scale structure formation as well. Moreover, the Planck mission will be able 
to measure the polarization of the CMB with unprecedented accuracy. I will show that polarization power 
spectra contain all information about the kinetics of hydrogen recombination and allow us to determine 
the parameters of the last scattering surface and the ionization history of the cosmic plasma at very high 
redshifts (z ~ 10 3 ). 

In the framework of the modern theory of the CMB anisotropy and polarization formation the kinetics 
of hydrogen recombination are assumed to be "standard" ones. The classical theory of the hydrogen re- 
combination was developed by Peebles (1968), Zel'dovich, Kurt and Sunuyaev (1968) for the pure baryonic 
cosmological model of the Universe and was generalized by Zabotin and Naselsky (1982), Jones and Wyse 
(1985), Seager, Sasselov and Scott (1999), Peebles, Seager and Hu (2000) for non-baryonic dark matter in 
the Universe. This standard model of recombination has been modified in various ways. 

Naselsky (1978) and Naselsky and Polnarev (1987) investigated the evaporation of primordial black holes 
as a possible source of the hydrogen recombination delay. Avelini et al. (2000), Battye et al. (2001), Landau 
et al. (2001) have pointed out that possible slow variations in time on the fundamental constants could be 
important for the ionization history of the cosmic plasma. Sarkar and Cooper (1983), Scott et al. (1991), 
Ellis et al. (1992), Adams et al. (1998) and Doroshkevich & Naselsky (2002) have investigated the influence 
of the decays of possible unstable particles on the kinetics of hydrogen recombination and have shown that 
these particles could be powerful sources for the distortions of the hydrogen ionization fraction. It is worth 
noting that all the models mentioned above show that the injection of non-thermal high energy electrons 
and/or photons produces the delay and distortions of the ionization history of the Universe starting from 
very high redshift z ~ 10 3 and down to redshift z ~ 5 — 10 when the galaxies started to form. These 
non-standard models on the ionization history of hydrogen recombination can be characterized by a few 
additional parameters which we needs to take into account in the framework of the best-fit cosmological 
parameter determination while using recent and future CMB anisotropy and polarization data. I call these 
parameters as "missing parameters" of the CMB anisotropy and polarization formation theory and I will 
show how important these parameters are for the CMB physics. The aim of my talk is to compare the 
possible manifestation of the more complicated ionization history of the Universe with the contemporary 
CMB anisotropy observational data (BOOMERANG, MAXIMA-1, CBI and VSA) in order to restrict some 
parameters of the models and to find the best-fit model of the hydrogen recombination, differing from the 
standard one. I'll show that a set of new parameters related with a more complicated ionization history of 
the Universe could play a significant role for the upcoming MAP and Planck data. 

2. Big Bang and the general properties of the CMB physics 

Our knowledge about the Universe is based on the experimental data and the well-developed theory, 
which we call the Bing Bang model: the inflation paradigm, baryogenesis, the theories of cosmological 
nucleosynthesis of the light chemical element, the formation of the CMB black body spectrum by Compton 
scattering of photons and electrons, the theories of hydrogen recombination and gravitational instability of 
the matter in an expanding Universe for the formation of large-scale structure of the Universe. These models 
serve as a background of the Big Bang theory and play an important role in the modern cosmology. Starting 
with a very hot and dense phase, the Universe expanded and cooled down. The black body spectra of the 
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Fig. 1. — The CMB image of the last scattering surface produced by the MAXIMA- 1 experiment. 



CMB and, in particular, the angular anisotropy and polarization preserve information about early epochs 
of cosmological evolution. When the temperature of the cosmic plasma dropped down to T ~ 3000/^ the 
free electrons recombined with protons into neutral hydrogen atoms and primordial hydrogen-helium plasma 
became transparent to the CMB photons. We call the corresponding redshift z. r ~ 10 3 as the redshift of the 
"last scattering surface", which is the epoch t r when the Thomson optical depth tt — 1. The recombination 
ceased when the optical depth dropped down to tt <C 1 during the corresponding time interval At r where 
At r /t r ~ |3/2Az r /z r | - 0.1. 

The questions arise as how exactly the cosmological plasma became neutral and how exactly it became 
transparent to the CMB photons. The answers to these questions lie in the physics of the CMB anisotropy 
and polarization formation. 

The angular anisotropy and polarization of the CMB radiation are related to primordial inhomogeneity of 
the cosmic plasma just at the epoch of the hydrogen recombination 1 Before this epoch (z > z r ) the optical 
depth of the plasma had been extremely high, so any ripples in the CMB temperature distribution caused by 
the primordial adiabatic perturbations of the metric, density and velocity of the plasma had been practically 
erased. After the epoch of recombination the mater component in the Universe became transparent and 
the CMB photons travel through the space and time from the last scattering surface to the observer. In 
Fig.l I show the image of the CMB anisotropy from a small path of the last scattering surface taken by the 
MAXIMA- 1 balloon-borne experiment. 

Actually, when we look at the last scattering surface, we see the processes of the structure formation just at 
the beginning. Small density and velocity perturbations are just traveling compressions and decompression 
of the photon-electrons gas similar to sound waves. These cosmological sound waves produce the ripples 
in the CMB map, which corresponds to the hot (white) and cold (black) spots in Fig.l. The cosmological 
sound waves are quite special. Before recombination the sound speed had been close to 1.7 x 10 4 km/s, but 
right after recombination it became 3 km/s! Such drastic transformation of the properties of the matter was 
described by Andre Sakharov in 1965 in connection with the matter density perturbations. The corresponding 
prediction by Sakharov is as follows. The amplitude of the perturbations after transition should preserve the 
modulation from the epoch before the transition should be preserved Because of the very high speed of sound 



see the lecture notes by A. Lasenby, J. Smoot and A. Melchiorri in this proceeding. 
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the typical wavelength of the perturbations is close to the size of the horizon at the epoch of recombination. 
It corresponds to the scale of 300 Mpc at the present moment of the expansion. The ratio between the 
spatial scale of the fundamental mode, which is <~ 300 Mpc, and the present horizon scale is exactly the 
angular measure of the CMB anisotropy perturbations, which is close to 30 arcmin. 

In discussion which follows, I will use a few characteristics of the CMB anisotropy. Firstly, let me introduce 
the fluctuations of the temperature AT in the spherical coordinate system with polar angle and azimuthal 
angle 0. 

We can decompose AT(0, <fi) as a function of and <j) with coefficients ai m on a sphere using spherical 
harmonics Yg m , 

AT(M)= ^ m Y lm {9,4>) (1) 

t m=-e 

The coefficients ae m are connected with the amplitudes and the phases of initial metric, density and ve- 
locity perturbations, which are the raw material for galaxy formation. This means that the initial conditions 
of the galaxies formation trough the coefficients ai m are related to the AT(0, 4>) distribution in the radio 
sky. 

Using aim we can define the power spectrum of the anisotropy and the perturbation of the temperature 
AT(£), which corresponds to the given value of the multiple number I. 

(a em a}, m ,) = Ce6e,t'S mt m' (2) 

As one can see from Eq (2), for random Gaussian processes there are no correlations between different 
harmonics and the phases of the multipoles are uniformly distributed in the range between and 2tt. 

In Figure 2 I show the numerical calculations of the CMB power spectrum for various cosmological models, 
which are reproduced from W. Hu's web site 2 . These two panels show the dependence of the CMB power 
spectrum on the baryonic density parameter fif,/i 2 and on the vacuum density parameter Q\h 2 . 

As one can see from Fig. 2 (a), the shape of the power spectrum looks like a set of peaks with fixed 

2 http: / /background. uchicago.edu/ whu/ 
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positions. For a flat universe the first peak corresponds to multiple number £ ~ 200. The second one 
£ ~ 400, the next one £ ~ 600 and so on. 

This is the manifestation of the sound waves in the power spectrum, which was predicted by Sakharov 
in 1965 and was widely discussed in the literature related to the formation of the CMB anisotropy and 
polarization power spectrum. In Fig. 2 (b) the dependence of the position and the amplitudes of the CMB 
anisotropy power spectrum peaks have shown for the A- and open CDM cosmological models for the given 
values of baryonic density Vt^h 2 = 0.0125. As one can see from this Figure, the CMB anisotropy power 
spectrum reflects directly all peculiarities of the ACDM cosmological model related with the geometry of the 
Universe at the redshifts z < 0.5. 

So, taking the above-mentioned properties into account, we can conclude that all peaks of the CMB 
anisotropy power spectrum form the basis for a new era of the so-called "precision cosmology" . Namely, major 
attentions of the COBE, BOOMERANG, MAXIMA-1, DASI, VSA, CBI experiments focus on determination 
of the power spectrum of the CMB anisotropy and testing of the statistical properties of the signals. In 
addition, all experimental teams have reported on the best-fitting cosmological models, which result in 
the most probable values of the baryonic matter density Cl),, the density of the cold dark matter £l c dm, 
the density of the vacuum Oa, possible tilt of the power spectrum of the adiabatic perturbation n s , the 
optical depth of reionization r r , the curvature parameter fix and so on, using most probable value of the 
Hubble constant h = 7?o/(100kmsec _1 Mpc _1 ) ~ 0.65 — 0.70. It is necessary to note that these values of 
cosmological parameters are obtained from very time-consuming methods, which include as a rule 6 or 11 
and more variables such as f2f>, f2 c , CIk . . . r r , ft and so on (e.g. see Tegmark et al. 2001). These parameters 
describe the most general properties of the Universe starting from very early epoch and down to the present 
moment of the expansion. 

An intriguing question is why we need to use 6, 11 or more parameters for determination for the most 
realistic cosmological model. The best explanation for high-dimensional space of the cosmological parameters 
has been suggested by Yu. N. Parijskij. Let us suppose that we have only one peak in the CMB power 
spectrum, i.e. the first peak at £ ~ 200 (see Fig. 2 a,b). For qualitative description of the peak we need 
three parameters: the amplitude of the peak, its width and location. For two peaks in the power spectrum, 
obviously, we need to introduce extra 3 parameters and for three peaks another 3 parameters and so on. From 
Fig. 2 we clearly see that at the multipole range £ < 1500 the number of peaks in the CMB power spectrum 
is N — 5 and thus corresponding number of independent parameters should be 15. So, generally speaking, 
the space of the cosmological parameters should be high-dimensional because of the peak-like shape of the 
CMB anisotropy power spectrum. 

Let me remind you very briefly of some general properties of the cosmological parameter determination 
from the CMB anisotropy data sets. The parameter dependence leads to the following definition of the CMB 
power spectrum for the initial adiabatic perturbations, 

C e = C (fit, O dm , O a , Ox , S1 p , fi™, n s , r, r r . . . [i]) (3) 

where Cldm is the dark matter density (including cold and warm components) scaled to the critical density of 
the matter p c , Q,k = 1 — J2j i s tne curvature parameter, which determines the flat, close or open model 
of the Universe, f2„ and fijj* are the density of massless and massive neutrino, n s is the spectral index for 
adiabatic perturbations (e.g. n s = 1 for the Harrison-Zeldovich power spectrum), n t is the corresponding 
spectral index for the gravitational waves, r — C*(£ — 2)/C s (£ — 2) is the ratio of the gravitational waves 
and scalar Cg power spectrum at £ — 2, t t is the optical depth of the hydrogen- helium reionization caused by 
the first quasar and galaxy activity. The notation [i] in Eq (3) is the main topic of my lecture, meaning that 
some additional parameters resulting from the history of hydrogen recombination could play a significant 
role in the transition of "the CMB power spectrum — > cosmological parameters" . Obviously, for the standard 
model of recombination we assume [i] = 0, which means that we do not need any additional parameters 
in Eq (3). If the kinetics of hydrogen recombination is distorted, however, by any sources of ionization or 
have some peculiarities of the spatial distribution of baryonic fraction, then it is necessary to include a few 
new parameters to the standard scheme. The question is how many parameters can characterize the more 



5 



complicated ionization history of hydrogen recombination at redshifts z < 1000 — 1500? In order to answer 
this question it is necessary to examine in details the "standard model" of hydrogen recombination. 

3. The standard recombination: an approximate model 

The recombination of free electrons with protons into neutral hydrogen atoms during the cosmological 
expansion means that atomic physics, and more specifically, the physics of hydrogen and helium atoms, 
contain some characteristics of the recombination process, which does not depend on the cosmology and can 
manifest themselves in a laboratory. Obviously, the corresponding atomic parameters are related with the 
structure of the electron energy levels in hydrogen atoms. 3 Fortunately for the physics of the CMB anisotropy 
formation, the physics of the electronic transition rates in hydrogen is well known from theoretical as well 
as experimental point of view. The frequency v n n i of a photon emitted or absorbed in a transition between 
two levels of energy in hydrogen with principle quantum numbers n, n is given by well-known equation from 
quantum mechanics 

E = hv = E n -E w , (4) 

where h is the Planck constant (not the dimensionlcss Hubble constant!), and E n = Ryn~ 2 , Ry = e 2 /2a = 
13,6eV, ao = h 2 /me 2 . To liberate a free electron from a neutral hydrogen atom one needs a photon at least 
the energy E — hv n n > — Xm where \n is the ionization potential from the initial state n. For n = 1 the 
corresponding potential is exactly I = Ry = 13,6eV and \n decreases as n~ 2 for n > 1. But it is exactly 
what we have in the cosmology just before the hydrogen recombination! 

As mentioned in the Introduction, the "hot" Big Bang model of the Universe indicates that the most 
widespread component in our Universe is the CMB photons, which has a perfect black-body spectrum and 
the corresponding number density of the photons n 7 ~ 420cm -3 at the present time of the cosmological 
evolution. In comparison, the number density of baryons n& — 2 x 10~ 7 (£!;,/i 2 /0.02) is significantly smaller 
than the CMB photons. The corresponding photon-baryon ratio s = n^/rib ~ 2 x 10 9 (£!h/i 2 /0.02) -1 is an 
excellent measure how hot the Universe is. Taking into account the temperature Tq = 2.728 from the COBE 
data for the black-body CMB spectrum (Bennet et al. 1996), we can find the number density of the Ly a 
quanta (I a — 10.2eV) in the Wien regime of the spectrum n(E) at the present time, 

/*°° ( hv \ 

n a = J n(E)dE ~ n 7 exp ^-yjrj (5) 

where I a — hv a is the energy of the Ly a quanta and k is the Boltzmann constant. If we reverse the 
cosmological expansion and trace back in time the temperature of the black-body radiation increases and in 
term of the redshifts it corresponds to T(z) = Tq(1 + z). When the redshift z is close to the redshift of the 
last scattering surface (z r ), then T(z r ) ~ 2700 K and the corresponding energy E = kT(z r ) ~ 0.3 eV. Thus, 
the number density of the Ly a quanta at z — z. r is n a (z r ) ~2x 10~ 15 n 7 and increase rapidly for the higher 
redshifts and reach the number density of baryons at z ~ 1500. For the Ly c quanta with E = I = 13.6eV 
the number density n c ~ n tt exp[-(/- I a )/kT(z)] <C n a at the redshifts z < 1500 and decreases rapidly 
when z — ► z r . So, taking properties of the electron transition in a hydrogen atom into account, one can 
conclude that photo-ionization of the neutral hydrogen atoms by the Ly a and Ly c quanta from the CMB 
power spectra should by a one of the main process of the recombination kinetics. 

One additional remark comes from the analysis of the photo-ionization cross-sections for resonant Ly a 
and Ly c quanta with neutral hydrogen atoms. In the framework of the theory of the CMB anisotropy 
and polarization formation the last scattering surface is characterized by the Thomson cross-section ctt = 
87r/3(e 2 /mc 2 ) 2 ~ 6.65 x 10~ 25 cm 2 , where e and m are the charge and the mass of an electron, and the 



3 I will discuss the kinetics of hydrogen recombination as a mean process of the formation of the last scattering surface. The 
kinetics of the cosmological helium recombination is included in the standard RECFAST code (Scager, Sassclov and Scott 1999) 
and I send off the reader to that paper. 
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corresponding optical depth tt — J aTn e cdt. For the photo-ionization of neutral hydrogen atoms which 
forms during the epoch of hydrogen recombination the corresponding optical depth 

a ph {u)n Hi cdt ~ 0ph - x~\t t \ z=Zt > r T \ z=Zr ~ 1 (6) 
a T 

where is the number densities of excited hydrogen atoms at different energy levels of the electrons, 
&ph{ v ) is the corresponding cross section of photo-ionization from the excited energy levels of the electrons 
to the continuum, c is the speed of light and x e ^ = n e /nji i . In the quantum theory of the spectral lines we 
have the following normalization of the photo-ionization cross sections (see for the reviews by Rybicki and 
Lightman 1979) 

2 

(J ph {v)dv = —fn^n' ( 7 ) 







where /„_>„' = f nn > is called the oscillator strength for the transition between states n and n . For example, 
the Ly a transition (n = l,n = 2) in hydrogen yields the /i^ value /i^ = 0.4162 (Rybicki and Lightman 
1979). Thus, we can conclude that for the Ly a transition the corresponding optical depth in the center of 
the line is approximately more that 10 orders of magnitude higher than the Thomson optical depth at the 
redshift z r . For the bound- free transition for hydrogen corresponding cross-section is given by the formula 
( Rybicki and Lightman, 1979): 



<7bf = 



where g is the Gaunt factor, a = e 2 /he = 1/137 is the fine-structure constant, a = h 2 /me 2 and hu) n = 
In, In = a 2 mc 2 /2n 2 is the ionization potential for the n — > continuum transition. For Ly c quanta the 
corresponding cross-section is 07 ~ oT 3 ot ~ 2 x 10 6 ctt and the optical depth is approximately 6 orders of 
magnitude higher than the Thomson optical depth at the moment of recombination. So, as one can see from 
the above- mentioned properties of the photo-ionization cross section, the CMB photons play an important 
role in the kinetics of hydrogen recombination in the Universe. Following Zcldovich, Kurt and Sunuyaev 
(1968) and Peebles (1968), let us discuss the quantitative model of this process during the epoch z < 1500. 

3.1. Approximate theory of hydrogen recombination 

As it was mentioned above in the expanding Universe the free electrons recombine with protons into 
neutral hydrogen atoms when the temperature of the matter-photon fluid drops down to T ~ 3000 — 4000 
K. As the first step of our investigation on the recombination process we have considered the distribution 
among the levels of a single atom in thermal equilibrium. We want to determine the distribution of an 
atomic species among its various stages of ionization, which is described by the well-known Saha equation. 
This equation describes the early stages of hydrogen recombination. Using the so-called three levels of the 
electron-energy model of hydrogen atom as shown in Fig. 3, we can start to describe the transition of the 
electron from the ground level to the continuum in thermal equilibrium (Rybicki and Lightman 1979). 

According to the Boltzmann law the differential number of ions at the ground level (N + (v)) with the free 
electrons at the velocity range between v and v + dv is 



N+(v) _ g 
No 9o 



exp 



I + \mv 

w 



2 1 



(9) 



where iVo and go are the number and statistical weight of the hydrogen atoms at the ground state, I = 13.6eU 
is the ionization potential, g = g$g e is the product of the statistical weight of the ion at the ground state 
(<7o") and the differential electron statistical weight g e — 2d 3 xd 3 p/h 3 . (Note that the factor 2 in g e reflects the 
two spin states of the electron.) The volume element d 3 x = n" 1 , where n e is the number density of the free 
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Fig. 3. — The three electron energy levels model of the hydrogen atom . 

electrons, and the volume element d 3 p is equivalent to 4wm 3 v 2 dv assuming isotropic velocity distribution. 
Using Eq. (9) after integration over all v we obtain 



where n p is the number density of protons and no is the number density of the hydrogen atoms at the ground 
level. As one can see from Eq. (10) the fraction of free electrons decreases very fast during the cosmological 
expansion because of the exponential factor (I/kT) » 1. However, I would like to put forward two important 
remarks while applying the Saha equation for the kinetics of cosmological hydrogen recombination. Firstly, 
it is necessary to note that during recombination the thermal equilibrium between black-body radiation 
and atoms has been affected by the deviation of the radiation spectrum from the Planckian form. Since 
recombination occurs for kT <C / the number density of photons with hv — I would be far smaller than 
the number density of electrons and protons. The range of the photon spectrum that is crucial for the 
recombination process actually changes substantially if recombining electrons proceed directly to the ground 
state. Secondly, it is necessary to take into account the cascading recombination process, which plays the most 
important role in the recombination kinetics and produce corresponding delay of hydrogen recombination 
relative to the "standard" fast Saha recombination kinetics. The importance of the cascading recombination 
process was pointed out by Zeldovich, Kurt and Sunuyaev (1968) and Peebles (1968). For the cascading 
recombination approximation we can neglect the direct recombination to the ground state but introduce the 
net recombination to the level (n, I), 



where a n i is the recombination coefficient to the (n, I) level of energy, n n i is the number density of the excited 
hydrogen atoms, (3 n i is the rate of photo-ionization of a hydrogen atom at this level and ni, is the number 
density of baryons at redshift z. Assuming that the process of recombination does not perturb significantly 
the black-body spectrum of the CMB at hv <C / range, the ratio of (3 n i to a ni is given by the Saha equation 
for n > 2, 




(10) 




(11) 




n>l 



(12) 
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where a c = X)n>i a ni an d In is the potential of ionization from the level n. In addition, it is necessary to 
take into account that in equilibrium the number density of atoms at the level n, I satisfies the following 



n n i = n 2s (2l + l)exp 



I a In 

kT 



(13) 



So, from Eq (1 1)-(13) we obtain the well-known relation for the number density of the free electrons (Peebles 
1968), 

d ( n„ 



- n b-n — 
at V ni, 



a c ni - l3 c G a n ls 



(14) 



where G a — n 2s /Tii s is the number of photons per mode in the Ly a resonance line. 

Let me describe briefly the evolution of the resonance Ly a line at the epoch of hydrogen recombination 
(Zcldovich, Kurt and Sunuyaev 1968, Peebles 1968). One of the most important results from the Ly a photo- 
ionization cross section is that the mean free time of the resonance recombination line of photons becomes 
very short, at the order ~ 30 seconds while the age of the Universe at the redshift z <~ 1500 is at the order of 
10 12 s. That means that the optical depth for the Ly a photons at the center of the line is extremely high, at 
the order of 10 10 and the radiation in the line is strongly coupled with the formed hydrogen atoms. On the 
other hand, the cosmological expansion of the Universe redshifted the Ly a photons from the center of the 
line and decrease rapidly the cross section of photo- ionization. To obtain the magnitude of the redshift effect, 
we can use the approximate equation for the spectrum of the electromagnetic radiation (Peebles 1968), 



d f TL 



(15) 



where 



n(y)dv 



(16) 



and H(t) is the Hubble parameter, i> m i n and ^ max are the minimal and maximal frequencies just below and 
above the center of the line, n(v) is the spectral number density of the photons in the line, R(t) is the net 
rate of production of the Ly a photons per unit volume during 2p — > Is transition. In order of magnitude, 
taking Eq. (16) into account, we have n a ~ n„ min Ai;, where Av is the Ly a line width and Avjv a ~ 10~ 5 . 
Thus, from Eq. (15) one can obtain (Peebles 1968) 

n Vrain = r^ max + RvaH^it) (17) 

Using dimensionless variable G — c 3 n{v) /Sirv 2 , we can change the last equation to the following form 

G a = G++R\lH- 1 {t)/8ir, (18) 

where G+ corresponds to n(j/ max ). Because of the intimate contact with the atoms the number of photons 
per mode satisfies the equation G + = exp[— (I a — I n )/kT] and for G a in Eq. (16) we get 



G 



cxp 



let In 

kT 



+ Rv a H- 1 (t). 



(19) 



The last point of the approximate model of hydrogen recombination is related to the estimation of the net 
recombination R(t), which defines the ionization fraction of the hydrogen. Thus 



R(t) = [a c n 2 e - p c G a n ls ] - A 2s ,u 



n 2s - n ls exp - 



I -la 
kT 



(20) 



where A 2s .is = 8.227sec 1 is the decay rate from 2s state. The excess of recombination over photo ionization 
results in the production of the Ly a quanta (the first term in Eq. (20)) or in a two-quantum decay of the 
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meta-stable 2s level of energy of the hydrogen atom (the last term in Eq. (20)). So, after simple algebra 
from Eq. (15)-(20) we obtain the well-known equation for the free electron number density (Peebles 1968): 



d_ 

dt 



C a c n\ - (3 c n u exp - 



I -la 

kT 



where 



C 



1 + KA2 S .isnis 



1 



(21) 



(22) 



K(A 2s ,i s + I3 c )n ls 

and K = X^ x H~ 1 (t)/8w. It is necessary to note that the approximate model of hydrogen recombination 
assumes that the temperature of the plasma T m during the recombination process is equivalent to black-body 
radiation temperature and the contribution of the helium atoms is negligible. A more general and detailed 
model of recombination was developed by Seager, Sasselov and Scott (1999), which includes the Compton 
interactions between electrons and photons and the corresponding deviation T m and T cm b, recombination of 
helium and ionization by the electrons. The qualitative theory of the hydrogen and helium recombination was 
summarize in the RECFAST program, which was incorporate in the CMBFAST code (Seljak and Zaldarriaga 
1996) for calculations of the power spectrum of the CMB anisotropy and polarization for various cosmological 
models. Below we will use the modification of the the RECFAST and CMBFAST codes to investigate more 
complicated than "standard" recombination regimes. 



4. Delayed and accelerated recombination 

As it was mentioned in Introduction and Section 3, the standard model of hydrogen recombination assumes 
that at the epoch z < 1500 the production of the resonant Ly a and Ly c photons is determined by the 
black-body spectra of the CMB and transitions of electrons into the neutral hydrogen atoms. Generally 
speaking, any decay of some relics of the cosmological evolution of matter in the Universe at the epoch 
of recombination are negligible. This assumption, however, needs to be proven by recent and future CMB 
experiments, especially by the Planck observation. 

As it was mentioned in Introduction a lot of non-standard models have been developed during last few 
years especially for the upcoming Planck experiment, which is characterized by unprecedent accuracy of 
the anisotropy and polarization power spectrum determination . The aim of this section is to compare the 
possible manifestation of the more complicated ionization history of the Universe with the contemporary 
CMB anisotropy observational data (BOOMERANG, MAXIMA- 1, CBI and VSA) in order to put constraints 
on some parameters of the models and to classify the models of hydrogen recombination, differing from the 
standard one. 

The basis for classification of any models of delayed and accelerated recombination is related with the 
model proposed by Peebles, Seager and Hu (2000), in which two independent parameters e a and £» are 
introduced as effectives of the non-thermal resonant Ly a and ionized-photon production before, during and 
after the period z ~ 10 3 , 

^=e a H(t)n H , (23) 

and 

^= £l H(t)n H , (24) 

where tih is the fraction of neutral hydrogen. In their model the parameters e a and Si are assumed to be 
constant in time, which means that corresponding number density of the Ly a and ionized-photon production 
per interval if -1 are An a = H~ 1 dn a /dt = e a nn and Am = H~ 1 dn i /dt = Sinn- As one can see, if e a <C 1 
and Ej«l then the corresponding distortions of the recombination kinetics should be small. It is natural 
to assume that the corresponding number density of free electrons produced by photo-ionization of the 
neutral hydrogen atoms is proportional to An, and An a and corresponds to the electron ionization fraction 
x e <~ Atii/riH ~ Ei or x e <~ An a /n H <~ e a for different stages of hydrogen recombination. However, the 
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concrete values and possible variation of the e a and Ei parameters in time depends on the energy injection 
spectra of the sources of the non-thermal resonant Ly a , ionized photons production and on the model of the 
source decay. 

Let us suppose that all resonant and ionized photons are the result of some massive unstable particles 
decay (Doroshkevich and Naselsky 2002), which include primordial black hole evaporation. For such particles 
the number density decreases in time as 

where nj is the number density of the X-particles, and Tx(nx,t) is the life time, which can be constant 
for the standard decay or as a function of cosmological time, if the production of the ionized and resonant 
quanta is related with annihilation of the x particles. After normalization nx = fix ( 2a ) , where ao is the 
scale factor of the Universe at some moment to, we can transform Eq. (25) to the following form 

^=e x (t)H(t)n H , (26) 

where 

ex = - (fe)- 1 exp f-'-^) ( (27) 

and fix, in is the initial number density of the A-particles at the moment to. 

I would like to point out that Eq. (26) has a general character. Obviously, we can include to the function 
ex(t) all models of the particle decay and then e a and Ei should be proportional to —ex(t)- As the result, 
they are functions of time (or redshift z). 

For illustration of the tight dependency e a , e, and £x(t) parameter, let me describe the model of the 
Super Heavy Dark Matter (SHDM) particle decay during the cosmological expansion. Such a model was 
discuss in connection with the origin of the Ultra High Energy Cosmic Rays (UHECRs) in our Galaxy in 
the framework of so-called Top-Down scenario. As was suggested by Berezinsky, Karshelriefi and Vilenkin 
(1997), Kuzmin and Rubakov (1998), Birkel and Sarkar (1998), the formation of such UHECRs can be 
related to decays of the various kinds of SHDM particles with masses Mx > 10 12 GeV. This model was 
widely discussed in connection with AGASA (Hayashida et al. 1994; Ave et al. 2000), Fly's Eye (Yoshida 
and Dai 1998; Abu-Zayyad et al. 2002) 4 , Havera Park (Lawrence, Reid and Watson 1991) experiments. 



4.1. Expected flux of high energy photons 

As is commonly believed, decays of SHDM particles into the high energy protons, photons, electron- 
positron pairs and neutrinos occurs through the production of quark-antiquark pairs (x — > q,q), which rapidly 
hadronize, generate two jets and transform the energy into hadrons (u>h ~5%) and pions (1 — ujh ^95%) 
(Blasi 1999). It can be expected that later on the energy is transformed mainly to high energy photons 
and neutrinos. This means that, for such decays of SHDM particles with 10 12 GeV< M x < 10 19 GeV, the 
UHECRs with energy E > 10 20 eV are dominated by photons and neutrinos (Blasi 1999). This conclusion 
can be tested with further observations of the UHECR fluxes at E > 10 20 eV. 

Probable energy losses of neutrinos are small (Bhattacharjee and Sigl 2000) , they are concentrated at 
the range E <~ Ei n j » Eqzk and could be responsible only for relatively rare observed events. But the 
comparison of expected and observed fluxes of high energy photons restricts some properties of SHDM 
particles and the Top-Down model as a whole ( Berezinsky, Blasi and Vilenkin 1998). At high redshifts the 
interaction of high energy photons with E ~ E in j 3> Eqzk with the CMB background (7uhecr + 7cmb — » 
e + + e~)lcads to formation of electromagnetic cascades (Greisen 1966; Zatsepin and Kuzmin 1966). 



4 see for details E. Loh's lecture in this proceeding 
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Fig. 4. — The left panel: the observed spectrum for the UHECR and Top-Down scenario predictions for 
different M x . The dash line is for M x = 10 n GeV, the dot line M x = 10 2 GeV, and the solid line M x = 
10 13 GeV. (From Sarkar and Toldra 2002. The right panel: the 7-rays spectra from UHECR in comparison 
with EGRET data. See for details in Sigl (2001). 



At small redshifts the efficiency of this interaction decreases and the evolution of such photons depends 
upon unknown factors such as the extragalactic magnetic field and properties of radio background. 

The expected domination of products of decay of the x-particle by the high energy pions and photons 
follows from quite general arguments. The cumulative contribution of the photons is bigger then the cu- 
mulative contribution of the hadrons by a factor of <~ oj^ 1 3> 1. This means that, at least for the galactic 
component, the cumulative observed flux of UHECRs at energy E > Eqzk is also dominated by galactic 
component of photons with E ~ Ei n j » Eqzk- The discrimination of high energy photon component of 
UHECRs can be considered as the crucial test for the Top-Down models. 

Summarizing available information about the photon component of UHECRs, Bhattacharjee and Sigl 
(2000) estimate the spectrum of injected photons for a decay of a single X-particle as follows: 

dN inj 0.6(2 -a) .U {2E i y a M x 

' E 7 < (28) 



dE 1 M x 0.9 \M X ) ' ' " 2 

where /„• is the fraction of the total energy of the jet carried by pions (total pion fraction in terms of of 
number of particles), < a. < 2 is the power index of the injected spectrum, M x is the mass of the SHDM 
particles. For the photon path length l y <~ 1 — lOMpc at E 1 ~ Eqzk in respect to the electron-positron pair 
creation on the extragalactic radio background (see Bhattacharjee and Sigl (2000); Protheroe and Bicrmann 
1996) we get for the photon flux jinj(E 7 ) at the observed energy E 1 : 

jinj ( El )^±l 7 (E 7 )n x ^ (29) 

where h x is the decay rate of the X particles. For the future calculation we will use the normalization of 
jinj (E 7 ) on the observable UHECR flux which corresponds to normalization on the decay rate h x at present 
time, t — t u (Bhattacharjee and Sigl 2000): 

h xfi 10- 46 cm- 3 s- 1 M 1 1 6 - a ex, (30) 

where 

^ 0.5 0.9 lOMpc / / E 2 obsJobs (E obs ) \ 
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and Mi 6 = M x /10 16 GeV, E le = £ ofcs /10 16 GeV, E obs and j obs (E obs ) are the observable range of the UHECR 
energy and flux, respectively. Note that normalization (30) does not depend on the nature of the x-particles. 

An additional possibility for non-equilibrium high energy photon injection to primordial hydrogen-helium 
plasma at the epoch of recombination comes from the theory of primordial black hole evaporation (Naselsky 
1978; Ivanov, Naselsky and Novikov 1994; Kotok and Naselsky 1998; Doroshkevich and Naselsky 2002) or 
relatively low massive particle decay Mx <C 10 12 GeV (Scott, Rees and Sciama 1991; Adams, Sarkar and 
Sciama 1998). For primordial black holes with masses M p bh ~ 10 13 gram the corresponding life time is 
Tpbh ~ t r and they evaporate just at the epoch of hydrogen recombination by emitting high energy particles 
with black body spectra with effective temperature T p bh ~ 10 13 — 10 14 K and energy -E p bh ~ fc^pbh ~ 1 — 10 
GeV (Hawking 1975). As one can see, this range of energy injection is sufficiently smaller than for the SHDM 
particles with Mx > 10 12 GeV. Note that for the models from Scott, Rees and Sciama (1991) and Adams, 
Sarkar and Sciama (1998) the characteristic energy of injected photons from neutrino decay is of the order 
10 — 100 eV, which is close to the Ly c — Ly a energy range. Taking the above-mentioned properties of the 
non-equilibrium photon injection into account it is natural to expect that practically all range of energy 
above Ly c — Ly a threshold should be very important for investigation of the non-standard recombination 
kinetics. 



4.2. Electromagnetic cascades at the epoch of the hydrogen recombination. 

The decays of X-particles cannot significantly distort the thermodynamic of the universe at high redshifts 
z > 10 3 but this injection of energy changes the kinetic of recombination at z <~ 10 3 that leads to observable 
distortions of the power spectra of CMB anisotropy and polarization. To evaluate these distortions we need 
firstly to follow the transformation of high energy injected particles to UV photons influenced directly the 
recombination process. 

The electromagnetic cascades are initiated by the ultra high energy jets and composed by photons, protons, 
electron- positrons and neutrino. At high redshifts, the cascades develops very rapidly via interaction 
with the CMB photons and pair creation (7uhecr + 7cmb — * e+ + e~), proton-photon pair production 
(Puhecr + 7cmb — > P + i + e+ + e~), inverse Compton scattering (cjjHECR + 7cmb — > e + 7 ), pair 
creation (guhecr + 7cmb — » e + e~ + e + + 7 ), and, for neutrino, electron- positron pair creation through 
the Z-resonance. As was shown by Berezinsky et al. (1990) and Protheroe et al. (1995), these processes 
result in the universal normalized spectrum of a cascade with a primary energy E 1 which can be written as 
follows: 

1 E a < E < E c ( 32 ) 
E r < E 



N 1 {E,E 1 ) = F{E,E 1 ) < 



F ( E > E -y) = o i w F /F ; > / EN 7 dE = E 1 



E-yE 2 

2 + ln(E c /E a ) 

where E c ~ 4.6 • 10 4 (1 + z)- 1 GeV, E a = 1.8 • 10 3 (1 + z)- 1 GeV. At the period of recombination z - 10 3 and 
for lesser redshifts both energies, E a and E c are larger then the limit of the electron-positron pair production 
E e + e - — 2m e — 1 MeV and the spectrum (32) describes both the energy distribution at E > E e + e - and 
the injection of UV photons with E <C E e + e - . However, the spectrum of these UV photons is distorted due 
to the interaction of photons with the hydrogen - helium plasma. 

In the range of less energy of photons, E < 2m e , and at higher redshifts, z > 10 4 , when equilibrium 
concentrations of Hi, He\ and Hen are small and their influence is negligible, the evolution of the spectrum 
of ultraviolet photons, N uv {E 1 z), occurs due to the injection of new UV photons and their redshift and 
Compton scattering. It is described by the transport equation (Berezinsky et al. 1990) 

dN uv 3N UV + _d_ ( Nuv ^] + Q(E,z) = Q 



dz l + z dE\ av dz) (1 + z)H 
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]_dE_ 
E~dz~ 



where the Hubble parameter is 



(l + z)H(z) 



E 



l + z 



Q{E,t)=h x N 1 {E,M x ) 



(33) 



H(z) = HoV^m(l + z) 3 + 1 - n ri 



(34) 

and n e oc (1 + z) 3 is the number density of electrons, so, /? 7 oc (1 + z) 3 l 2 E. Here Q(E, t) is considered as an 
external source of UV radiation. 

The general solution of equation (33) is (Doroshkevich and Naselsky 2002): 



N uv (z) 



E(x) 



l + z 



E(z) 1 



Q{x) E 2 {x) 
H{x) E 2 (z) 

- 2 ME,z) 



1 



1 + x 



l + z 



dx 
1 + x' 

5/2 

- 1 



(35) 



where the maximal redshift, z mxi in (35) is defined by the condition E(x) = 2m e c 2 . 
As is seen from (35), the Compton scattering dominates for E » 30keV, when 



p 7 (E,z)=U 



E 



0.3 hn h 



1 



m e c 2 V 0.02 



3/2 



10 3 

VTTz~e T (z 



H(z)p 7 (E,z) ~ E 5 / 2 
For the most interesting energy range, E <C 30keV, (3j(E, z) <C 1, we get again 

2«x(z) 



N uv {E(z), z) 



3 ff(z) 



■N 1 {E,E 1 =M x /2). 



(36) 
(37) 

(38) 



It is necessary to note that the UV photons energy spectra Eq. (38) reflect directly the simple idea that 
practically all energy of the SHDM relics transforms to the energy of the 7 quanta at the range E < 100 
GeV (at z — 0)(see Fig4.b) and than the soft part of the spectra forms due to electromagnetic cascades. 
At higher redshift, especially at z <~ 10 3 the electromagnetic cascades from the 7-range of the spectra are 
very effective for the UV photon production because of the high energy density of the CMB (12 orders of 
magnitude higher than the present energy density) and because of increasing of the mean energy of the CMB 
photons , which is in order magnitude equivalent to 1 eV. 

4.3. Delay of the hydrogen recombination due to SHDM decay 

Let me briefly discuss the relationship between initial spectra of the SHDM decay and corresponding 
parameters e a and Si of the hydrogen excitation and ionization . As it was mentioned in the previous 
subsection, the spectra of UV photons can be normalized to the energy density of injection Q(t) from the 
SHDM decay as 



N{E) ~ n- 



Q 



-3/2 



(39) 



E 

'E 2 a [2 + \n(E c /E a )} \E~ a 

where 77 is the effectives of the energy transformation from the E ~ range to the 7 range of the 

photons energy spectra. Let us normalize the energy density of injection Q(t) (eVcm~ 3 s -1 ) to the present 
day value Q{t u ) 

Q(t) = w(t)Q(t u ) (40) 
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where t u is the age of the Universe at z = and w(t) oc ex according to Eq. (25). Integration Eq. (39) over all 
energies of the photons corresponds to the net of the energy density injection e 7 = / dEEN(E) = u)(t)Q(t u )- 
But for the energy range E > I a the net of the ionized UV photons is given by integral J dEN(E): 

,41, 



E a {z) 

Using normalization for the energy density e^H^ 1 ~ eegret to the EGRET energy density cegret — 
410~ 7 cVcm~ 3 (Sigl 2001) of the 7 quanta at the present moment of the cosmological expansion, we can 
estimate the function u(t). Namely, assuming 



we can find that 



w(t) oct 4 -P oc (l + z)i( 4 - p ) (42) 

* * *« * «*) (y^f^y) (1 + z)*™ oc (1 + (43) 

where = (//^(z)) 1 / 2 ~3x 10~ 6 \/TTI. The power index p = 1 correspond the model involving release 
of x-particles from topological defects, such as ordinary cosmic strings, necklaces and magnetic monopoles 
(see for review by Sigl 2001). The model p = 2 corresponds to the SHDM exponential decay (see Eq. (25)) 
with t x 3> tjj- From Eq.(43) one can see that for p = 4/3 the e% and e a parameters does not depend on the 
rcdshift, which is exactly the case of the model by Peebles, Seager and Hu (2000). 

One additional possibility comes from the model with tx < t u . This model could be related with 
primordial black holes evaporation or all SHDM relics decay, for which the corresponding life time is less 
than the present age of the Universe. Note that for tx — 0.1t u (and corresponding rcdshift zx ~ 6 — 7 ) 
the SHDM relic decay before or just at the galaxy formation epoch and the observable flux of the UHECRs 
mostly should relate with extragalactic background of the particles, which is characterized by very small 
anisotropies. For such a model we need to modify Eq. (43) as £ mo d{z) — £(z)@(z) for z >• Zx, where 

9(zx)=exp[(l + zx) i ] , (44) 

and zx corresponds to the tx- 

In Fig. 5 I show the dependence of the parameter on the redshift z for various models of the SHDM 
decay (for different p and Q(zx))- 

As one can see from Fig. 5, some of the SHDM relics can distort the kinetics of hydrogen recombination 
quite significantly {e ai £i > 10~ 3 — 10~ 2 for z < 1500), for some of the SHDM relics the perturbation of 
the ionization history of the Universe is small and practically unobservable. However, all energy injection 
during the epoch of hydrogen recombination would by partially absorbed by the hydrogen atoms and produce 
corresponding delay of recombination or ionization of the neutrals. 

The question is we can have the acceleration of the hydrogen recombination and what we need to know 
about the properties of the cosmological plasma beyond the standard model? 



4.4. Acceleration of recombination 

In the framework of non-standard recombination models, it was shown that all peculiarities of the hydrogen 
photo-ionization kinetics can be described in terms of e a , Si parameters, which correspond to additional non- 
equilibrium Lyaipha and Ly c photon injection. In the previous section it was shown that for different models of 
the SHDM particles decay both parameters e a , e» are positive and hence the delay of hydrogen recombination 
occurs. In terms of the e a , £j parameters any acceleration of recombination means that one or both of them 
are negative. However, any decay of the SHDM relics do not produce "negative" energy injection and, taking 
into account the black-body CMB photons, acceleration of the recombination means that we can remove 
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Fig. 5. — The dependence of the e» parameters versus z for different models of the SHDM decay. Solid 
line corresponds to p = 0,Q(zx) = 1 model, dot line is p = l,Q(zx) = 1 model, dash dot line is the 
model p = 2, @(zx) = 1 and dash line and dash dabble dot line are the models p = 2, Q(zx) = 10 6 and 
p = 2, Q(zx) = 3 x 10 5 correspondingly. 



some of the photons from interaction with the formed neutral hydrogen atoms, which is impossible. Docs 
it mean that for given values of the cosmological parameters (il b ,Q m . . .etc.), which determines the general 
properties of the model of the Universe, we can only have the delay of hydrogen recombination? To answer 
this question, we need to focus on Eq. (21), in which the net of free electron number density decreases during 
recombination which is proportional to the electron and proton number density. The non-linear dependence 
of the ionization fraction on baryonic number density is the key to the construction of the accelerated model 
of recombination! 

The basic idea of accelerated recombination is very transparent (Naselsky and Novikov 2002). Suppose 
that, instead of pure adiabatic model of initial metric, density and velocity perturbation of matter we have the 
mix between adiabatic and isocurvature (or isotcmpcraturc) perturbation. If at the very small mass scales of 
isocurvature perturbations, e.g. those comparable with 10 2 — 1O 5 M0, the distribution of the baryons is very 
non-uniform like clouds, beyond which scale perturbations follow the standard adiabatic properties, then 
the hydrogen recombination inside the clouds goes faster than outside the clouds. This indicates that for 
adiabatic perturbations which produce anisotropies and polarization of the CMB, the hydrogen ionization 
differs from the standard one, which corresponds to the mean value of the baryonic matter density. Non- 
trivial is that such a modification of the ionization history of the hydrogen in the cloudy baryonic model can 
be described in terms of e a < and £, > parameters as for the delayed recombination model (Naselsky 
and Novikov 2002). All of these delayed and accelerated recombination models could be very important for 
reconstruction of the best-fitting cosmological model parameters from the modern CMB observational data 
and can distort the optimal values. 

5. Non-standard recombination and its manifestation in the CMB anisotropy and polarization 
power spectrum 

For illustration of the introduced above classification of the hydrogen recombination regimes we will use 
two cosmological models with the following parameters. The Model 1: Q h h 2 — 0.022, f2 c = 0.125, = 
0, fl\ = 0.7, T r = 0.1, h = 0.7, n s = 1, s a = £; = 0. This model has the best-fitting of the CBI observational 
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data (Mason et al 2002). The Model 2: tt b h 2 = 0.032, fl c = 0.115, fl K = 0, Q, A = 0.7, r r = 0.1, ft = 0.7, n s = 
1, £ a = £i = 0. This model differs from the previous one on the baryonic and CDM fraction densities and it is 
non-optimal according to the CBI data. Using RECFAST code we can find the differences in the ionization 
fraction of hydrogen plotted in Fig. 6. 

Fig. 7 is the plot of the ionization fraction of the hydrogen in the clumpy baryonic model (Model 3), 
which corresponds to the parameters of the Model 1, but instead of O^ft 2 = 0.022 we have to use the mean 
value of the baryon density (ilfc)ft 2 = 0.022, assuming that density contrast between inner and outer regions 
£ = pin/ Pout = 11, the fraction of volume in clouds / = 0.1 and the mass fraction Z — £//[l + (£— 1)/] = 0.5. 

The mass spectrum of a clouds means to be close to S(M — M c i), where 10 2 M Q < M c i < 10 6 M Q . For 
this model the mean ionization fraction is related with ionization fractions inside Xi n and outside x out clouds 
as imcan = Xi n Z + x ou t(l — Z) (Naselsky and Novikov 2002). In Fig. 8 we show the difference between £ mean 
and x\ normalized to (x\ + x moan )/2. 

As one can see from Fig. 8, at the range 700 < z < 1500 we have acceleration of recombination while at 
the range z < 700 we have delay of recombination. The last model, Model 4, which we would like to include 
in our analysis is the Model 2 plus additional sources of delay and ionization, corresponding to the following 
values of the e a ,Si parameters: 



£ a (x) — ax exp 



1 + z 
l + z x 



1 + z 
1 + z x 



£i{z) = -s a (z) (45) 
a 

where a ~ 0.3, (3 ~ 0.04 and z x ~ 10 3 . As one can see from Eq(6) at z ~ 10 3 the variation of the e ai £i~ 
parameters are very slow and e a ~ a, Si ~ /3for comparison with the parameters by Peebles, Seager and Hu 
(2000). In Fig. 9 I plot the ionization fraction in Model 4 in comparison with ionization fraction in the Model 
1 and the Model 2. As one can see the delay of recombination in the Model 4 produce 20 — 30% deviation 
of the ionization fraction at the redshift z ~ 10 3 and very significant deviations at the redshift z < 10 3 . 

It is worth noting that all the Models 1-4 clearly illustrate the different delayed and accelerated regimes, 
discussed at the beginning of the this section. The question is how sensitive the CMB anisotropy and 
polarization power spectrum on different models of ionization history of hydrogen and baryonic density? 
Below I will show that combination Q(,,£ a ,£j is very important for estimation of the correct value of the 
baryonic fraction density from the CMB data. 

5.1. Anisotropy and polarization as a test on ionization regimes 

In order to compare the CMB anisotropy and polarization power spectrum in the Models 1-4 we have 
to use modification of the CMBFAST code (Seljak and Zaldarriaga 1996), by taking into account more 
complicated ionization history of the plasma. All the models which include the late reionization correspond 
to the standard CMBFAST option for a given value of the optical depth r r . This value of r r is related with 
the ionization fraction x = 1 at the redshift of reionization z r as 

2 /3 /O u2\ -2/3 /r) ,2x1/3 



where Y p is the helium mass fraction of the matter. In Fig. 10 we plot the CMB power spectrum for the 
Models 1-4 in comparison with the BOOMERANG, MAXIMA-1, and CBI data at the multipole range 
I < 2000 in which the possible manifestation of the Sunuyaev-Zel'dovich effect is not important. 

As one can see from Fig. 10 all the models look very similar to each other. However, in the Table 1 we show 
the values of the x 2 -parameter for all the models. From the Table 1 we seen that the Model 3 with baryonic 
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Fig. 6. — The rate of ionization for the Model 1 and 2. Solid line corresponds to the Model 1 (xi), dash line 
Model 2 (x 2 ), dash-dot line the ratio (x\ — x<2)jx\ . The reionization epoch are not included. 
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Fig. 7. — The rate of ionization for the clumpy baryonic model. Dash-dot line corresponds to the ionization 
fraction for the inner regions at. Dot line is the ionization fraction for the outer regions. Solid line the mean 
fraction of the ionization. Dash line is the ionization fraction for the model with ilf, = 0.045. 



Table 1: \ 2 f° r observations and models 



obs/mod 


BOOM 


MAXIMA 


CBIM1 


CBIM2 


VSA 


Model 1 


11.2 


14.1 


3.20 


7.15 


7.08 


Model 2 


34.6 


17.7 


2.75 


7.51 


13.0 


Model 3 


8.69 


11.3 


1.60 


5.62 


6.59 


Model 4 


27.7 


15.7 


1.87 


5.15 


10.0 
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Fig. 8. — The difference of the rates of ionization between the clumpy baryonic cosmological model and 
without clumps. 




Fig. 9. — The rates of ionization and the differences of the rates of ionization between the Model 1, 2 and 
4. Solid line is the ionization fraction in the Model 1. Dotted and short-dash line are the Model 2 and 4, 
respectively, which look very similar. Long-dash line corresponds to the ratio 2(x4 — x\)/(x\ +x&). Dash-dot 
line corresponds to 2(x± — x 2 )/(x2 + X4). 
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clouds has an excellent agreement with the CBI data while the Model 1 is characterized by the value of the 
Xi = 2x1- Note that the Model 1 is the best-fitting model for the CBI data without any assumption on 
more complicated ionization history of cosmological hydrogen recombination. The Model 4 corresponds to 
the ilbh 2 parameter 1.5 times bigger than the Model 1. However, this model also has an excellent agreement 
with all the CBI observational data. 

From the practical point of view, the Models 1 , 3 and 4 consistent with the CBI data and more complicated 
ionization history of the hydrogen recombination manifest themselves as new sources of degeneracy of the 
cosmological parameters, namely, the baryonic density of matter. This conclusion is very important for 
the upcoming MAP and Planck data. In order to characterize the differences between the models and to 
compare their with the sensitivity of the upcoming Planck experiment we plot in Fig. 11 the functions 

D h3 (£) = 2 [d(£) - C 3 (£)} I \C x {i) + C 3 (£)} ; 

D 2A M = 2 [C 2 (£) - C 2 (£)] I [C 2 {I) + C 2 {£)] ; (47) 

for the multipole range 2 < I < 2500. Obviously, these functions D lt3 {£) and D 2 ^{£) is necessary to compare 
with the errors of the Ci extraction for the Planck mission. 

As one can see from Fig. 11, the range of the Di^^) and D 2 ^(£) values is close to the 5% . Note that 
it can easily drop down to 1 — 2% using smaller values of the parameters of the non-standard ionization 
history e a ,Ei and corresponding parameters £, Z, zx, m and so on. These parameters play a role of the 
"missing parameters" in the model of the CMB anisotropy formation and significantly increase the number 
of the standard parameters fib, fl\,- ■ ■ £Ir- It is well known that for analysis of the degeneracy of some of the 
cosmological parameters the polarization measurements are very useful. In Fig. 7 we plot the polarization 
power spectrum AT 2 (£) = £(£+l)/2wC p (£) for the Models 1-4. As one can see from this Figure, the difference 
of the ATp(£) is very small at the range £ > 200 and quite significant for the Model 1 and 2 relatively the 
Models 3 and 4 because of the bump at £ <~ 5 — 10, which is related with reionization of the hydrogen at 
the epoch of formation of galaxies and clusters of galaxies. Some useful information can be obtained at the 
range 20 < £ < 200 in which the manifestation of the Model 4 is significant. This Model shows that the 
optical depth r,. = 0.1 can be obtained without distortions of the low multipole part of the power spectrum 
at £ <~ 5 — 10. However, for the Planck mission one of the most important source of uncertainties at the 
multipole range I < 200 is the the cosmic variance. The corresponding errors of the AT 2 (£) estimation at this 

range is 5C P {£)/C P {£) ~ (/ sk y^) _1/2 ^ 9 (£/200) _1/2 (/ sky /0.65) _1/2 %. 5 For / sky /0.65 ~ 1 and £ ~ 2000 
the cosmic variance limit on 5C(£)/C(£) is 3% at 68% confidence level, which means that all the peculiarities 
of the anisotropy power spectra below this limit should be unobservablc. For the CMB polarization, as it 
follows from Fig. 8 all features of the power spectrum can be observed by the Planck and probably, the 
low multipole part of the power spectrum can be observed by the MAP mission. So future experiments for 
detection of polarization can be very useful for investigation of the ionization history of the cosmic plasma 
and possible distortions of the kinetics of hydrogen recombination by the different sources. 

6. Conclusion 

The precision of the measurements on the CMB anisotropy and polarization, especially being expected 
from the MAP and Planck missions, allows us to distinguish the cosmological models under consideration 
and to discuss 10 — 20% distortions of the standard recombination process occurring at redshiftz ~ 10 3 . 

In my talk I have compared available data from observations of the CMB anisotropies with several cosmo- 
logical models with different kinetics of hydrogen recombination, which is caused by possible external sources 
of the energy injection and possible small-scale clustering of the baryonic component in the clouds. It was 
shown that such a sources can produce the delay or acceleration of hydrogen recombination and increases 
the number of parameters used to fit observed power spectra of the CMB anisotropy and polarization. 



We assume that instrumental noise and systematics should be close to the cosmic variance limit . 
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Fig. 10. — The CMB power spectrum for the Models 1-4. Solid line corresponds to the Model 1, dash line 
the Model 2, dash-dot line the Model 3, and long-dash line the Model 4. 
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Fig. 11. — The differences in the CMB power spectrum for the Models 1-4. Solid line corresponds to the 
Di^(£) function and dash line I?2,4(^)- Thick solid line is the estimated errors for the Planck . 
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Fig. 12. — The CMB polarization power spectrum for the Models 1-4. Solid line corresponds to the Model 
1, dot line the Model 2, dash line the Model 3, and dash-dot line the Model 4. 
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Fig. 13. — The differences of the CMB polarization power spectrum for the Models 1-3 and the Models 2-4. 
Definition for residuals are similar to the definition for the anisotropy. Solid line is for D\^{t) and dash line 
D2,a{(-), but for polarization. Thick solid line is the estimated errors for the Planck . 
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The most interesting conclusion coming from the investigation of the non-standard models of hydrogen 
recombination is that the CMB polarization measurements will play an important role for the estimation of 
the "missing parameters" from the upcoming Planck data. 
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